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© Echo canceller. 



© An echo canceller modem having a fast training echo canceller is described. The echo cancellation 
coefficients are computed by taking correlations between a complex, periodic training signal sequence and a 
real component of the corresponding echo signal. The modem includes a receiver circuit which detects a signal 
on a channel possibly including an echo. An echo canceller estimates the real component of the echo signal. 
Training circuitry applies the complex training sequence to the channel and takes correlations between the 
training sequence and the real component of the corresponding echo signal. In the preferred arrangement, the 
modem also includes a computational element for computing the period of the periodic sequence, generating the 
complex periodic training sequence in real time and for computing a phase roll frequency based upon the 
computed echo cancellation coefficients. 
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This invention relates to data communication equipment or modems provided with echo cancellation 
features. 

Data communication equipment (DCEs), or modems, are devices used to transmit and receive binary 
data over a communication channel. One category of DCEs, sometimes referred to as full-duplex modems, 

5 is capable of performing the functions of transmitting and receiving simultaneously. When the modem is 
transmitting and receiving simultaneously over a two-wire communication link (e.g., transmission over the 
switched telephone network), an echo of the transmitted signal is commonly present in the signal received 
from the remote modem. If the transmitted and received signals occupy the same frequency band, it is 
necessary to remove the echo signal, in order to reliably detect the data sent by the remote modem. 

10 The echo signal typically has a near echo and a far echo component. The near echo is generated by 
the imperfect hybrid couplers in the local modem and the near-end telephone central office. The far echo, 
on the other hand, is mainly generated by the hybrid couplers in the remote central office and the remote 
modem. The far echo is delayed in time relative to the near echo. When this delay can be substantial, the 
echo canceller is often broken into a near-echo and a far-echo canceller component which are also 

75 separated by a delay. 

High-speed modems typically use bandwidth-efficient modulation schemes such as quadrature modula- 
tion. In such systems, the binary data is first mapped into a sequence of complex signal points (symbols) 
chosen from a constellation with a finite number of points. The real-valued transmitted signal carries 
information about this complex sequence. 
20 Two-wire, full-duplex, high-speed modems, such as the standard V.32 voiceband modem specified by 
the CCITT, come equipped with adaptive echo cancellers which are capable of nearly eliminating the 
echoes of the transmitted signal. An echo canceller is typically implemented as a transversal filter which 
consists of a tapped-delay line, and a series of variable complex-valued tap coefficients. The inputs to the 
tapped-delay line are the aforementioned complex signal points. These are appropriately weighted by the 
25 tap coefficients to generate as output the real part of the weighted running sum. This represents an 
approximation of the received real-valued echo signal. The echoes are cancelled by subtracting this 
estimated echo signal from the real received signal. 

Echo cancellers which are implemented as a transversal filter with a complex-input and a real output 
are often referred to as Nyquist echo cancellers. Nyquist echo cancellers often consist of a near canceller 
30 and a far echo canceller. One realization of Nyquist echo cancellers is described by S. Weinstein in the 
U.S. Patent 4,1 31 ,767 (reissue Re31 ,253). 

An echo canceller is typically trained in the absence of the remote signal, during an initialization or 
training period which occurs prior to data transmission. In many echo cancellers, the transversal filter is 
trained using the least mean-square (LMS) algorithm. In an LMS algorithm the tap coefficients are 
35 continually adjusted to remove any correlation between the complex input symbols and the residual 
received signal which remains after echo cancellation. However, the time required to accurately train an 
echo canceller in this manner can be very long, particularly in modems which employ echo cancellers with 
long transversal filters. 

In the past, fast training methods have been discussed for echo cancellers whose input and output are 
40 either both real- or both complex-valued. One method using an ordinary periodic chirp sequence was 
disclosed by T. Kamitake in IEEE Proc. of ICC'84 (pp. 360-364, May 1984, Amsterdam, Holland) in a paper 
entitled "Fast Start-up of an Echo Canceller in a 2-wire Full-duplex Modem". A similar method using a 
pseudo-random shift-register sequence was later described by V. Kanchan and E. Gibson in IEEE Trans, on 
ASSP (Vol. ASSP-86, No. 7, pp. 1008-1010, Jul. 1988) in a paper entitled "Measurement of Echo Path 
45 Response". These methods are not applicable to Nyquist echo cancellers with a complex input and a real 
output. 

In IEEE Proc. of GLOBECOM'87 (pp. 1950-1954, Nov. 1987, Tokyo, Japan), J. M. Cioffi proposed a 
method for fast simultaneous training of both the near and the far echo cancellers of a Nyquist echo 
canceller in a paper entitled "A Fast Echo Canceller Initialization Method for the CCITT V.32 Modem." This 

50 method is based on the discrete Fourier transform (DFT) and uses a real periodic pseudo-noise sequence, 
however, it assumes a perfect Hilbert filter (transformer) in the transmitter which is generally not realizable. ' 

When the echo path can be substantially modelled as a linear filter, a transversal filter can effectively 
reconstruct and substantially cancel the echo signal. However, in certain instances, the far echo may also 
contain a small amount of frequency offset, also referred to as phase roll, which complicates echo 

55 cancellation. Some more sophisticated echo cancellers also include phase-roll compensation circuitry which 
can track the phase variations in the far echo and thereby remove its detrimental effect. Typically, the 
phase-roll compensation circuitry includes a phase-locked loop (PLL) to acquire the phase-roll frequency 
and phase during the training period. 

2 
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Summary of the Invention 

One aspect of the invention is defined by the features of claim 1 . 

A second aspect of the invention is defined by the features of claim 7. 
5 Further embodiments are defined in the dependent claims. 

Practical embodiments of our modem directly estimate echo canceller coefficients without having to use 
a few thousand bauds of training time to iteratively approach the optimum values. Furthermore, we have 
found the methods now taught for computing the echo canceller coefficients to be accurate and computa- 
tionally efficient as compared to other methods of computing such coefficients. In addition, the period of 
io time during which echo response data must be gathered for the computations can generally be reduced to 
a period which is significantly less than the total echo delay of the channel. Consequently, practical 
embodiments of our modem yield very fast echo canceller training, significantly faster than is achievable in 
previous modems. 

We have found that practical embodiments of our modems offer the additional advantage of producing 
75 echo canceller coefficients which are very close to their optimal values and which may result in a residual 
echo which is within 2 db of its optimal value. Thus, further training using alternative methods such as a 
least-mean-squares (LMS) algorithm may not be necessary. 

We can estimate the complex transfer function of a linear system when only the real output of the 
system is available. 

20 Other advantages and features will become apparent from the following description of preferred 
embodiments, with reference to the accompanying drawings, in which:- 
Fig. 1 is a block diagram of an echo cancellation modem; 

Fig. 2 depicts the autocorrelation properties of a complex periodic training sequence which is used to 
train the echo canceller shown in Fig. 1 ; 
25 Fig. 3 depicts the power spectrum of the complex periodic sequence which has the autocorrelation 
properties shown in Fig. 2; 

Fig. 4 is a flow chart of an algorithm for determining the period of the complex periodic training 
sequence; 

Figs. 5a-c depict the echo response of a representative channel over which the modem shown in Fig. 1 
30 communicates; 

Fig. 6 is a flow chart of another algorithm for determining the period of the complex periodic training 
sequence; 

Fig. 7 is a block diagram of a portion of the echo cancellation modem depicted in Fig. 1, depicting 
details of the trainer; 

35 Fig. 8 is a flow chart of an algorithm for determining the phase roll frequency and the power of the far 
echo; and 

Fig. 9 is a high level block diagram of a modem which embodies the invention. 
Referring to Fig. 1, in an echo cancellation modem 2, which may be of a Nyquist echo canceller type, a 
scrambler/encoder 6 receives a data bit stream over an input line 4. The scrambler/encoder 6 randomizes 

40 the data bits to ensure that any bit pattern is as likely to occur as any other bit pattern and it encodes the 
bits according to the coding system being used to generate a first sequence of complex symbols. A 
modulator 8, using a carrier signal 10 of frequency f c , modulates the output of the scrambler/encoder 6 to 
generate a digital transmit signal 12 consisting of a second sequence of complex symbols. Next, a 
transmitter filter 14, a digital-to-analog converter 16 and a low pass filter 18 convert the digital transmit 

45 signal to an analog signal 20 which is ready for transmission over a channel 22 to a remote device, not 
shown. A hybrid coupler 24 couples the analog signal 20 to the channel 22. 

The hybrid coupler 24 also accepts signals sent over the channel 22 to the modem 2 by the remote 
device and sends them to a receiver portion of the modem 2 as a received signal 26. During full-duplex 
communication over the channel 22, the received signal 26 includes both a near echo and a far echo, which 

so may have an associated phase roll. A band pass filter 28 processes the received signal 26 and an analog- 
to-digital converter 30 converts it to a real digital received signal 32. The analog-to-digital converter 30 
produces samples at times kT/M, where T is the baud interval of the local modem's transmitter, M is an 
integer chosen such that no aliasing will occur in the received signal after sampling, and k is a sampling 
interval index. (In the embodiment described herein, M is chosen to be equal to 3.) A combiner 34 

55 combines the samples of the real digital received signal 32 with corresponding samples of an estimated 
echo signal 36 to generate an echo-cancelled signal 38. Finally, a receiver 40 processes the echo-cancelled 
signal 38 to produce a received data bit stream 42 corresponding to the bit stream sent by the remote 
device. 

3 
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An echo canceller module 44 within the modem 2 generates the estimated echo signal 36. In the echo 
canceller module 44, a delay line 46 receives the complex transmit signal 12 and produces two groups of 
signals 48a-b which are separated in time from each other by a bulk delay. The group 48a corresponds to 
the near echo signal and includes a plurality of complex symbols each of which are delayed in time relative 

5 to each other so that they span the near echo signal; whereas the group 48b corresponds to the far echo 
signal and includes a plurality of other complex symbols each of which are delayed in time relative to each 
other so that they span the far echo signal. The group of delayed complex symbols 48a, corresponding to 
shorter periods of delay, are processed by a near echo canceller 50 and the group of delayed complex 
symbols 48b, corresponding to longer periods of delay, are processed by a far echo canceller 52. Both the 

10 near echo canceller 50 and the far echo canceller 52 include variable coefficients which are trained so that 
the echo cancellers 50 and 52 produce outputs which are accurate estimates of the near echo and the far 
echo, respectively. A summer 54 combines the outputs of the near echo canceller 50 and the far echo 
canceller 52 to generate the estimated echo signal 36. 

The echo canceller module 44 also includes a phase roll compensator (PRC) 56 which controls the far 

75 echo canceller 52. After being trained, the PRC 56 introduces a phase roll into the output of the far echo 
canceller 52 which approximates the phase roll of the far echo in the received signal 26. 

Training of the variable coefficients of the echo canceller module 44 is controlled by a training signal 
generator 58 and a trainer 60. During training, which may occur prior to data communication or during a 
period interrupting data communication, the training signal generator 58 generates a special echo training 

20 sequence that is sent to the remote device and the trainer 60 monitors the resulting real digital received 
signal 32. The remote device remains silent during this period of echo training, thus, the resulting real 
digital received signal 32 basically consists only of near and far echoes. Based upon the monitored echo 
signal, the trainer 60 computes echo canceller coefficients which set the variable coefficients in the near 
echo canceller 50, the far echo canceller 52 and the PRC 56 so that the echo canceller 44 generates an 

25 estimated echo signal 36 closely approximating the monitored echo. The echo canceller module 44 uses 
the echo-cancelled signal 38 as an error signal indicating how accurately the echo canceller module 44 has 
estimated the echo using the coefficients estimated by the trainer 60. Further fine adjustment of the echo 
canceller coefficients may be incorporated using the echo-cancelled signal 38 and the well known least- 
mean-square (LMS) adaptive algorithm. 

30 The training generator 58 generates a complex periodic sequence, p(n) = p r (n)+jpj(n), having a period of 
2L and the following autocorrelation properties: 



35 (n) =0,L 

Mod 2L 



40 



45 



50 



k=0 " lv) 



elsewhere 



2L^l J Mod 2L 

P i ( n )® 2L Pi(")= £Pi< k > x Pi«" + k) Mod2L )H-A < n > Mod2L = L 
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10 elsewhere 
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PrC")®^ p.(n) = Zp i (k)xp i ((n + k) Mod2L )=0 for all n. CV > 

k=o 



where: 

k is an integer delay index; 
n is the time index; 
55 L is the one half of the period of the periodic sequence; 

® 2L denotes the circular correlation operation of length 2L; 
A is a positive real number; and 

()mod2l means the modulo 2L operation required by the circular correlation. 
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As indicated the real and imaginary parts of the training sequence are orthogonal, i.e., the cross- 
correlation of length 2L between p r (n) and pKn) is zero for all n. The properties of the autocorrelation of the 
real and imaginary parts of p(n) are illustrated in Fig. 2. And the corresponding power spectrum of such a 
training sequence has the form illustrated in Fig. 3. Note that the power spectrum of p r (n) and p { (n) are the 
5 2L-point discrete Fourier transforms of the correlations given by Eqs. (1 ) and (2), above. 

Using a training sequence having the above characteristics, the trainer 60 calculates the echo canceller 
coefficients for the near echo and the far echo cancellers 50 and 52 by computing correlations between a 
period of the complex training sequence 12 and the real digital received signal 32. To appreciate that the 
computed correlations yield the desired values for the echo canceller coefficients, it is useful to review the 
w mathematics describing the operation of the system. 

The echo channel can be considered as a linear system with a complex impulse response c(t). If the 
transmitted data symbols are denoted as d(nT), where T is the baud time interval and n is the baud interval 
index, then the output of the system sampled at every T, y(nT), is related to the input, d(nT) as follows: 

75 y(nT) = Re[c(nT)od(nT)] (4) 

where © denotes the linear convolution operation. In this example, y(nT) represents the real sampled digital 
received signal 32. Actually, if the received signal 26 is sampled at a rate 3/T, as it is in the above- 
described embodiment, then the output signal y(t) may be written as y(nT + mT/3), where m = 0, 1 .2. This 
20 can be more simply denoted as y m (n). Using the same convention, c(nT + mT/3) can be rewritten as c m (n) 
which is the impulse response of the echo channel sampled at nT + mT/3 and equation (4) becomes: 

Ym(n) = Re[c m (n)0d(n)] (5) 

25 Substituting the complex training sequence p(n) for d(n) results in the following: 

y m (n) = Re[c m (n)op(n)] (6) 

After the channel is fully excited by the training sequence p(n), the real echo y m (n) is also periodic with a 
30 period of 2L, if the channel is time-invariant, and then the correlation between the real echo and the 
sequence is: 

y m (n)®2LP(n) = Re[c m (n)op(n)]® 2U p(n) (7) 
35 Separating c m (n) and p(n) into their real and imaginary parts yields: 
ym(n)® 2L p(n) = [c mr (n)op r (n) - c ml (n)©p,(n)]® 2 L[p r (n) + jp { (n)] (8) 

where c mr (n) and c mj (n) are the real and imaginary parts of c m (n), respectively, and j is V-1. Since the real 
40 and imaginary parts of the chosen p(n) are orthogonal under the circular correlation of length 2L, the above 
equation reduces to: 

y m (n)®2LP(n) = c mr (n)o[p r (n)® 2 i_Pr(n)] - jc mi (n)0[pj(n)0 2 LPi(n)] (9) 
45 y m (n)®2Lp(n) = c mr ((n) MO D 2l) + c mr ((n-L) M0 D 2i_) - jc mj ((n) MO D 2l) + jc mi ((n-L) MO D 2l) (10) 
ym(n)® 2L p(n) = c* m ((n) MOD 2l) + c m ((n-L) MOD 2l) (11) 

where * denotes the complex conjugate operation. Because of the periodicity, only the correlation for 
so 0£n£2L-1 need be considered. From Eq. (11) it is apparent that the correlation results for 0^n^2L-1 contain 

two duplicated (the first of which is complex conjugated) versions of c m (n), which have a span no greater 

than L. Hence, one of the versions provides enough information. 

In general terms, the properties of the training sequence, as shown in Figs. 2 and 3 and as specified in 

Eqs. 1-3, have the following significance. Note that the objective is to determine the impulse response of 
55 the echo channel which, of course, theoretically can be directly determined by stimulating the echo channel 

with an impulse. There is at least one problem, however, with using an impulse. Its energy is too 

concentrated in time and the impulse peak may drive the signal into nonlinear regions of the channel where 

the data signal does not normally operate during normal communications. By distributing the power evenly 



5 



BNSDOCID: <EP. 



06301 19A2_I_> 



EP 0 630 119 A2 



over time, as is implied by Fig. 3, one can transmit more energy with the training sequence and thus can 
more fully excite the channel without driving the transmitted signal into its nonlinear regions of operation. 
Moreover, by using a training sequence whose autocorrelation function is an impulse train, as shown in Fig. 
2, one can still readily and directly determine the impulse response of the echo channel by simply 
5 correlating the training sequence with the echo response corresponding to the training sequence. That is, 
the resulting function, namely, y m (n)® 2 LP(n), corresponds to the echo channel response that one would 
obtain by stimulating the channel with the impulse train represented by the autocorrelation of p(n). The 
orthogonality property of the real and imaginary parts of the training sequence simply guarantees that the 
real and imaginary parts do not interfere with each other, when the real and imaginary parts of the 

io coefficients are estimated. 

According to the correlation properties given by Eqs. (1) and (2), every half period of the correlation in 
Eq. (9) yields an estimate of the sampled echo channel response c m (n), if the total span of the echo is less 
than the half period L. Thus, it is verified that by using the periodic sequence characterized by Eqs. (1), (2) 
and (3), the real and imaginary parts of c m (n) can be obtained by simply performing circular correlation 

75 between a 2L-long segment of the samples of the real digital received signal 32 and the real and imaginary 
parts of one period of the periodic sequence p(n). There will be a sign difference between the imaginary 
parts of the estimated and the actual responses, if the correlation from the first half period is used. Since 
the real digital received signal is periodic, an alternative way to estimate the echo canceller coefficients c m - 
(n) is to perform a linear correlation of length 2L between the real and imaginary parts of one period of the 

20 periodic sequence p(n) and a 3L-long segment of the real received signal 32. 

For a Nyquist echo canceller with three subcancellers, the correlation processing is repeated for each 
incoming T/3 sample. Namely, Eq. (9) is repeated for m=0,1,2. Thus, there are three incoming samples per 
baud and three subcanceller coefficients are obtained per baud interval. 

A number of periodic sequences have the properties described in Eqs. (1) through (3). However, it is 

25 desirable that the sequence have a small peak-to-RMS (root-mean-square) ratio so that there is a lower risk 
of driving the transmitted signal into regions of nonlinear operation. A periodic sequence which has this 
desirable property is one for which the phase of its discrete fourier transform obeys a square law of 
frequency. Such a sequence is defined as follows: 
(1) For 0^n<L, and if 

30 (a) L is even, then 



35 




P r (n) = 



40 



P,(n) = 




V2L 



45 



(b) L is odd, then 



50 



55 
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L-1 

2 2 

,kn k 



1 + 2 f c °sI 2 ,(if + J5 rT ) ] 



L-1 
2 



p i( n) = 



75 (2) For L^n<2L, then the second half of the sequence is the complex conjugate of the first half of the 
sequence, namely: 

Pr(n) = Pr(n-L), and (16) 

20 Pi (n) = - Pi (n-L) (17) 

Another sequence which has a slightly higher peak-to-RMS ratio but is easier to compute in real time 
since fewer cosine values are required is: 
(1)For 0£n<L, 

25 

1 + cos((n+— )k)+2 ^ cos[ — ^ L ) 

so Pr( n) = — : 

INT ^ L 2 



35 



cos[(n + ±±±)k) * 2 g cos(2n[^^ip^ + 
P i (")= '-^ 



40 (2) For L^n<2L, 

p r (n) = Pr(n-L), and (20) 
Pi(n) = -Pi(n-L) (21) 

45 

where INT[(L-1)/2] is the largest integer which is less than or equal to [(L-1)/2]. 
Modem 2 implements the algorithm illustrated in Fig. 4 to determine the value of L. The algorithm 
determines the period which must be used in constructing the periodic sequence so that the training 
sequence produces interleaved, non-overlapping near and far echoes when the channel is fully excited. 
50 What this means can be more clearly understood with the aid of Figs. 5a-c. 

Fig. 5a shows the echo response of the channel to a single impulse at t = 0. The near echo occurs first, 
followed by the far echo delayed in time by B. Note that the span of the near echo and span of the far echo 
is less than or equal to Ni and N 2 , respectively. When the channel is excited by the impulse train shown in 
Fig. 5b, the echo response is as shown in Fig. 5c. Until the far echo appears, the echo response consists of 
55 a train of near echoes, one occurring every L baud. After B bauds have elapsed, the channel becomes fully 
excited and the far echoes begin to appear, also separated by a period of L. To utilize the near and far 
echoes thus generated for computing the echo canceller coefficients, it is desirable that the near and far 
echoes are non-overlapping as shown in Fig. 5c. In addition, it is desirable that the total time period, L, over 
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which the two echoes occur should be as short as possible so that the computations involve the minimum 
number of bauds. The algorithm shown in Fig. 4 selects the L which satisfies these two criteria. 

During the period of echo canceller training but prior to transmitting the special echo training sequence, 
the modem measures the far echo round trip delay, B, introduced by the channel 22 (step 100, Fig. 4). The 

5 delay, B, is equal to the length of time required for the far echo to come back after the training signal is 
transmitted. One such method for determining the far echo delay is described in the CCITT V.32 standard. 

Note that Ni and N 2 are usually determined at the time of designing the modem. Generally, they are 
empirically determined by studying the typical characteristics of channels. It is preferable to select Ni and 
N 2 so that they are upper bounds on the duration of the corresponding echoes which one expects to 

10 receive. 

After modem 2 has measured the far echo delay, B, it computes a variable, k, which is defined as 
follows: 

k = INT[(B - Ni)/(Ni + N 2 )] (step 110) 

75 

where INT[x] means the largest integer that is less than or equal to x. 

Next, modem 2 tests k to determine whether it is greater than zero (step 1 1 5). If k is greater than zero, 
the modem computes two additional variables, pi and p 2 , defined as follows: 

20 pi = (B + N 2 )/(k + 1) (step 120) 

p 2 = (B - Ni)/k (step 130) 

Then, in step 140, modem 2 determines whether an integer exists in the range R defined as follows: 

25 

pi ^ R ^ p 2 

If there is an integer within the range R, modem 2 sets L equal to the smallest integer within the range 
(step 150). That is, if pi is an integer, then L = pi ; otherwise, L = INT[pi +1]. The modem then branches to 
30 step 160 where it computes f, defined as follows: 

f = B-kL 

The variable f indicates the location of the far echo within the period L. The far echo is located between f 
35 and f + N 2 -1. Using this information, modem 2 assigns the correlation calculations for each baud in the 
interval of length L to the corresponding component of the echo. 

In step 140, if there is no integer within the range R, then modem 2 branches to step 170 where it 
decrements k by one and then returns to step 120. Steps 120 through 140 are repeated until the range R 
includes an integer. 

40 If in step 115 k is not greater than zero, then modem 2 tests if k equals zero (step 116). If k equals 
zero, modem 2 sets L equal to B + N 2 and then branches to step 160 to compute f. However, if k is less 
than zero, then modem 2 sets L equal to Ni +N 2l sets f equal to Ni (step 118) and then exits the algorithm 
(step 119). 

The algorithm yields an L which lies within the following interval: 

45 

Ni + N 2 ^ L ^ 2(Ni + N 2 ) - 1, 

for any practical values of Ni , N 2 , and B. 

For finite precision implementations of the algorithm, the selected value for L may be incorrect or non- 
50 optimal. If L is incorrect, the near and far echoes may overlap. Whereas, if L is non-optimal, the obtained 
period is not the smallest possible period. These errors are due to round-off error effects. The critical 
stages of the algorithm are in step 140, where the range R is tested for the presence of integers, and in 
step 110, where the initial value of k is computed. To avoid the problems associated with round-off errors, 
the algorithm may be modifed as shown in Fig. 6. 
55 Steps 200 and 215 correspond to steps 100 and 115, respectively, of Fig. 4. Step 210 is slightly 
different from step 110, namely, k = INT[(B-Ni )/(Ni +N 2 ) + So] where 5 0 is a small positive correction 
number. In step 215 if k is greater than zero, modem 2 computes pi (step 220) and p 2 (step 230) which are 
defined differently from what was described above. Namely, 

8 
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pi = (B + N 2 )/(k + 1) + 5i (step 220) 
p 2 = (B - Ni )/k + 5i (step 230) 

5 

where 5i is a small correction number. 

Then, in step 240, modem 2 determines whether an integer exists in the range R as previously defined, 

i.e.: 

10 pi £ R £ P2 

For this test, however, pi is treated as an integer if its fractional part is less than a threshold 5 2 . 

As before, if there is an integer within the range R, modem 2 sets L equal to the smallest integer within 
the range (step 250). That is, if pi is an integer, then L = pi; otherwise, L = !NT[pi] + 1. The modem then 
15 branches to step 260 where it computes f, as previously defined, i.e.: 

f = B-kL 

Finally, modem 2 determines if Ni^f<L-N 2 (step 280). If it is, the algorithm stops (step 219). 

20 If an integer does not exist within the range (step 240) or if f does not lie within the defined interval 
(step 280), modem 2 branches to step 270 where it decrements k by one. Then, it branches back to step 
215 to repeat the steps until an L which satifies the conditions of the algorithm is found. 

In step 215, if k is not greater than zero, modem 2 branches to step 216 which corresponds to step 116 
of the algorithm shown in Fig. 4. Indeed, the sequence followed after branching to step 216 is the same as 

25 that described for Fig. 4. That is, modem 2 implements steps 217 and 218 which are the same as steps 
117 and 118, respectively, in Fig. 4. 

The correction numbers, 5 0 and 5i , and the threshold, 5 2 , depend upon the word length and the 
rounding scheme used in the computations. It is preferable that the word length used in computation should 
be long enough. For example, when 16 bit fixed-point arithmetic is used, double precision computations are 

30 preferable. In this case, experiments have determined that when 5o =2~ 15 , 61 =2~ u and 5 2 =2~ 12 , then the 
algorithm computes a correct L and step 280 may not be necessary. 

The initial values of the coefficients of the near and far echo cancellers are computed as follows. The 
required half-period L of the special training sequence is computed based on the value of B as described 
above. The trainer 60 then computes, in real time, the required periodic complex training sequence 

35 according to the value of L. The computed training complex sequence is then transmitted by the transmitter 
for at least k f + 4 half-periods, or (k' + 4)L bauds, where k* = k if kZO or k* = 0 otherwise, and k is obtained 
during the determination of L. After at least k f + 1 half-periods, or (k* + 1)L bauds, of the training sequence 
have been transmitted, the trainer 60 stores the real received digital signal samples 32 in 3 buffers 70, 72 
and 74, shown in Fig. 7, which are T-spaced delay lines that are each at least 2L samples long. Another 2L 

40 bauds later, after each buffer has received 2L samples, and while the transmitter continues sending the 
periodic sequence, trainer 60, using correlators 76a-f, starts to perform correlations of the input samples 
with the real and imaginary parts of one period of the periodic training sequence. The correlations can 
either be circular correlations if the same samples in the 2L long buffers 70, 72 and 74 are being used, or 
they can be linear correlations if first-in-first-out (FIFO) buffers that receive 3 new samples per baud are 

45 employed instead. The real coefficients are directly obtained from the correlation results, while the 
imaginary part of the coefficients are equal to the correlation results multiplied by (-1) k Assuming three 
complex coefficients, one for each subcanceller, are computed per baud, a time period of Ni bauds will be 
used to compute all the coefficients of the near echo canceller 50 (see Fig. 1). The coefficients of the far 
echo canceller 52 are computed in the same way as computing the near canceller coefficients. The 

so coefficients of the far canceller 52 can be computed in N 2 bauds. 

The modem uses the calculated far echo canceller coefficients to set the variable coefficients of PRC 
56. PRC 56 may employ digital phase-locked loop (PLL) technology such as that described in U.S. Patent 
Number 4,813,073 entitled "Echo Cancellation", issued on March 14, 1989, the disclosure of which is to be 
regarded as incorporated herein by reference. To train the PLL in PRC 56, it is necessary to accurately 

55 estimate the phase-roll frequency and to measure the far echo power. Modem 2 implements the algorithm 
illustrated in Fig. 8 to compute both of these characteristics. 

Assuming that the frequency of the far echo phase-roll « p radians/s, the far echo is modulated by exp- 
[jopt]. The far echo received at time Ti will have a phase-rotation <*3 p Ty radians relative to the far echo 
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without phase-roll. If the frequency *> p is relatively low, such that in a 2L long segment, all the received 
signal samples have approximately the same phase rotation, the correlation property between the training 
sequence and the received signal described above will approximately still hold. An estimated far canceller 
coefficient using the received signal at the time Ti will be equal to the coefficient estimated without phase- 

5 roll multiplied by expU«pTi ]. Similarly, the same far echo canceller coefficient estimated at a later time T 2 is 
multiplied by expU^ p T 2 ]. Thus, the phase-roll frequency <* p can be estimated by dividing the angle 
difference, * = (T 2 -Ti )o> p , between these two far echo canceller coefficients by the time difference T2-T1. 
The estimation procedure is described below. 

First, modem 2 computes at least two sets of far echo canceller coefficients using the correlation 

io technique described above (step 300). A first set of coefficients is computed corresponding to time Ti and 
a second set of coefficients is computed corresponding to time T 2 , which occurs D baud after time Ti . In 
practice, it is convenient to choose D = L, although D can be another positive integer value. Then, in step 
310, modem 2 computes the estimated phase difference between the two sets of calculated far echo 
canceller coefficients. 

75 One way of determining the estimated phase difference is by computing the sine of the phase 
difference, </> m>n , between each of the corresponding coefficients in the two sets of calculated far echo 
canceller coefficients. 

Sin 4> m , n is related to the calculated coefficients in the following way: 



20 



25 



30 



45 



Im[c m (n+D) ]Re[c m (n) J-Im[c m (n) ]Re[c ffl (n+D) ] 
sin <f> m n = — — — — — 

I <=„(") | |c B (n+D) | 

This is approximated as follows: 



c ffi i (n+D) c m r (n) -c m * (n) c m r (n+D) 



where c mir (n) and c m>j (n) are the real and imaginary parts of c m (n), m=0,1,2, respectively, |c m (n)| is the 
35 magnitude of c m (n) and it is assumed that the magnitudes of the two calculated coefficients are equal. 
If the absolute value of the angle £ m , n is small, e.g. less than 20 degrees, then 

4>m,n * sin 

40 Theoretically, the angle between any pair of far echo coefficients should be the same. However, noise 
and other interferences may undercut the validity of this relationship. Therefore, to reduce the effect of any 
noise or interference and improve the accuracy of the estimate of the phase roll frequency, the estimated 
angle may be averaged over all of the pairs of far echo coefficients. That is, in step 310 modem 2 may 
compute the following: 



AVG«>) = E n>m a m . n * m , n 



where E n>m a m>n = 1, in order to assure that AVG(<£) is an unbiased estimate of the angle. 

To obtain an optimal weight a m , n - it is preferable to set a m , n equal to c 2 mi «(n) + c 2 mj (n)/Z min [c? mti (n) + c 2 m>r 
50 (n)]. When this is done, note that AVG(<£) equals: 



^m.nCc^ (n) c mi (n+D) -c mi (n) (n+D) ] 



55 ^ n [c 2 n(i (n)+c 2 B(r (n)] 



10 



BNSDOCID: <EP 06301 1 9A2J_> 



EP0 630 119 A2 



It may be preferable to estimate the average angle by computing the just-cited equation rather than by 
computing the angle for each coefficient and then averaging all of the computed angles. The reason is that 
the former approach involves only one division; whereas the later approach involves many divisions. Since 
division is a time consuming and inefficient process when using commercially available digital signal 
5 processors, it is desirable to minimize the number of divisions. 

After the estimated phase difference is computed, modem 2 computes an estimate of the average 
phase roll frequency by dividing AVG(<f>) by DT = T 2 -Ti (step 320). 

Finally, in step 330 modem 2 determines the power of the far echo by computing the following 
relationship: 

10 

P f = E |Cf (i)d(n-i) = |E f (n)| 2 = |d(n-i)| 2 ri|c f (i)| 2 

where the overbar denotes the ensemble average operation and d(n) is the data symbol which is 
transmitted. The value of P| is then used for optimal scaling of the PLL coefficients in PRC 56. 

75 The modem may be implemented by a multiple processor architecture, as shown in Fig. 9. That is, it 
has a general host processor 62, which performs overall control and data movement functions; a signal 
processing element 64, which performs the functions of the transmitter and echo canceller, including the 
implementation of the algorithms described above; and another signal processing element 66, which 
performs the functions of the receiver 40. A modem generally of this type is described in our European 

20 Patent Specification Publication No. EP-A-01 54551 , the disclosure of which is to be regarded as incor- 
porated herein by reference. 

Claims 

25 1. A modem for transmitting to and receiving signals from a remote device over a channel (22), the 
received signal including a near echo and a far echo separated by a delay; characterised in that the 
modem comprises: an echo canceller (44); a training generator (60) adapted for operatively generating 
a periodic training sequence with a symbol time interval T to train the echo canceller, the periodic 
sequence having a period equal to an integer multiple of a variable L, the near echo having a span of 

30 about NiT and the far echo having a span of about N2T, N1 and N2 being integers; and the near echo 
and the far echo being separated by a delay, BT; and a computational element adapted for operatively 
determining the value of the variable L by: selecting an integer K which is no greater than (B-N1)/- 
(N1 +N 2 ), and setting L equal to an integer existing in an interval R which is substantially defined as 
follows: [(B + N 2 )/(K + 1 )] f£ R <£ [(B-N1 )/K]. 

35 

2. A modem according to Claim 1 , further characterised in that the integer K is the largest integer that is 
no greater than (B-N1 )/(Ni + N 2 ). 

3. A modem according to Claim 1or 2 further characterised in that L is selected to be the smallest integer 
40 existing in the interval R. 

4. A modem according to any of Claims 1 to 3 further characterised in that the period is equal to 2L. 

5. A modem according to any preceding claim, wherein the far echo may have a phase roll, further 
45 characterised in that the echo canceller module is adapted to be controlled by specifying an estimate 

of the phase roll frequency and by using the said computed correlations to determine a set of variable 
coefficients in the echo canceller, and in that the modem further comprises: a computational element 
adapted for operatively: computing a difference in phase between a first set and a second set of echo 
canceller coefficients, both sets having been generated by the trainer module, the first set of echo 
50 canceller coefficients corresponding to a time T1 and the second set of echo canceller coefficients 

corresponding to a later time T 2 ; and 

dividing the computed phase difference by the time difference, T 2 -Ti, to arrive at the estimate of the 
phase roll frequency. 

55 6. A modem according to claim 5, further characterised in that said computational element is adapted to 
compute the computed phase difference as a weighted average of the phase differences between 
corresponding elements of the first and second sets of echo canceller coefficients. 
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7. A method for computing a period used to construct a periodic training sequence for training an echo 
canceller, the training sequence being applied to a channel and resulting in a corresponding echo 
signal including a near echo having a span of about Ni T and a far echo having a span of about N 2 T, 
Ni and N 2 being integers and T being the symbol time interval of the training sequence, and the near 
echo and the far echo being separated by a delay , BT; the method being characterised in comprising 
the steps of: 

selecting an integer K which is no greater than (B-Ni)/(N! + N 2 ), and setting the period equal to an 
integer multiple of an integer L existing in an interval R which is substantially defined as follows: 

[(B + N 2 )/(K + 1 )] <> R £ [(B- Ni )/K] 

a A method according to Claim 7, further characterised in that the integer K is the largest integer that is 
no greater than (B-Ni )/(Ni + N 2 ). 

9. A method according to Claim 7 or 8 further characterised in that L is selected to be the smallest integer 
existing in the interval R. 

10. A method according to Claim 7, 8 or 9 further characterised in that the period is equal to 2L. 

11. A method where according to Claim 7, 8, 9 or 10 wherein the far echo may have a phase roll, further 
characterised in that the operation of the echo canceller is controlled by specifying an estimate of the 
phase roll frequency and by using the said computed correlations to determine the echo canceller 
coefficients, and in that the method further comprises the steps of: 

computing the difference in phase between a first set and a second set of echo canceller coefficients, 
both sets having been generated by the trainer module, the first set of echo canceller coefficients 
corresponding to a time Ti and the second set of echo canceller coefficients corresponding to a later 
time T 2 ; and 

dividing the computed phase difference by the time difference, T 2 - Ti , to arrive at the estimate of the 
phase roll frequency. 

12. A method to Claim 11, further characterised in that the computed phase difference is a weighted 
average of the phase differences between corresponding elements of the first and second sets of the 
echo canceller coefficients. 
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© An echo canceller modem having a fast training 
echo canceller is described. The echo cancellation 
coefficients are computed by taking correlations be- 
tween a complex, periodic training signal sequence 
and a real component of the corresponding echo 
signal. The modem includes a receiver circuit which 
detects a signal on a channel possibly including an . 
echo. An echo canceller estimates the real compo- 
nent of the echo signal. Training circuitry applies the 
complex training sequence to the channel and takes 
correlations between the training sequence and the 
real component of the corresponding echo signal. In 
the preferred arrangement, the modem also includes 
a computational element for computing the period of 
the periodic sequence, generating the complex peri- 
odic training sequence in real time and for comput- 
ing a phase roll frequency based upon the computed 
echo cancellation coefficients. 
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